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Continuous gravitational-wave (CW) signals such as emitted by spinning neutron stars are an important tar- 
get class for current detectors. However, the enormous computational demand prohibits fully coherent broad- 
band all-sky searches for prior unknown CW sources over wide ranges of parameter space and for yearlong 
observation times. More efficient hierarchical "semicoherent" search strategies divide the data into segments 
much shorter than one year, which are analyzed coherently; then detection statistics from different segments 
are combined incoherently. To optimally perform the incoherent combination, understanding of the underlying 
parameter-space structure is requisite. This problem is addressed here by using new coordinates on the param- 
eter space, which yield the first analytical parameter-space metric for the incoherent combination step. This 
semicoherent metric applies to broadband all-sky surveys (also embedding directed searches at fixed sky posi- 
tion) for isolated CW sources. Furthermore, the additional metric resolution attained through the combination 
of segments is studied. From the search parameters (sky position, frequency, and frequency derivatives), solely 
the metric resolution in the frequency derivatives is found to significantly increase with the number of segments. 

PACS numbers: 04.80.Nn, 95.55.Ym, 95.75.-z, 97.60.Gb 



I. INTRODUCTION 

Direct detection of gravitational waves would not only 
validate Einstein's theory of general relativity but also con- 
stitute an important new astronomical tool. Continuous 
gravitational-wave (CW) signals are expected, for instance, 
from rapidly rotating neutron stars through various emission 
mechanisms ||T|-2]- Most such stars are anticipated to be 
electromagnetically invisible, but might be detected and stud- 
ied via gravitational waves. With current Earth-based detec- 
tors, such as LIGO [61, numerous efforts are underway to 
search for CW sources [7. -101 . and observational upper limits 
have already allowed one to constrain the physics of neutron 
stars imfTJll. 

The expected CW signals are extremely weak. Hence de- 
tection requires very sensitive data analysis techniques to ex- 
tract these signals from detector noise. In the work of [TTl 
a powerful method has been derived which is based on the 
principle of maximum likelihood detection, leading to coher- 
ent matched filtering. The CW signals are quasimonochro- 
matic with a slowly changing intrinsic frequency. For a ter- 
restrial detector, the Earth's motion relative to the solar sys- 
tem barycenter (SSB) generates a Doppler modulation in am- 
plitude and phase of the waveform. As shown in lfT3l . the 
parameters describing the signal's amplitude variation may be 
analytically eliminated by maximizing the coherent matched- 
filtering statistic. Thus, one only searches over the remain- 
ing parameters describing the signal's phase: the source's 
sky location, frequency, and frequency derivatives ("spin- 
downs"). The so-obtained coherent detection statistic is called 
the J^-statistic, which can also be generalized to employ mul- 
tiple detector data streams lfT4l . 

Finite computational resources are what imposes severe 



limits on the sensitivity of broadband all-sky searches for 
prior unknown CW sources |fT3] \\5 \. In the fully coherent 
J^-statistic approach, one must convolve the full data set with 
signal waveforms (templates) corresponding to all possible 
sources. But the number of templates required to cover the 
search parameter space increases as a high power of the co- 
herent integration time. For year-long data sets, searching a 
realistic range of parameter space would demand more com- 
puting power than available on Earth ifTSIfTSl . Therefore, a 
fully coherent search is restricted to much shorter integration 
times. 

A more efficient analysis of data sets which span yearlong 
time intervals is achieved via less expensive hierarchical semi- 
coherent methods lfT6l - [T9l . In such a method, the data are di- 
vided into segments of duration T, where T is much smaller 
than one year. This allows one to use a coarse grid of tem- 
plates, on which the J^-statistic is calculated for each segment. 
Then the J^-statistics from all segments (or statistics derived 
from F) are incoherently combined using a common ^ne grid 
of templates. Since phase information is discarded between 
segments, this latter step is called incoherent and thus the 
search methodology as a whole is refeiTed to as semicoher- 
ent. 

An important long-standing problem in these semicoher- 
ent methods has been the design of, and link between, the 
coarse and fine grids. To address this problem it is essential to 
understand the underlying parameter-space structure. In this 
context, the geometric approach has been proven to be very 
useful, introducing a metric on parameter space, as first done 
in 1201 I2TI . The key quantity in this respect is called mis- 
match M., which is the fractional loss in expected J^-statistic 
(or sum of expected J^-statistics in the incoherent step) for a 
given signal p^j at a nearby grid point p in phase-parameter 
space V. Taylor-expanding M (to quadratic order) at p^j in 
the differences of p^- and p defines the parameter-space met- 
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While the parameter-space metric for coherent searches has 



been examined in detail [ 15 , 22 , 23], the metric for semicoher- 
ent searches is comparably much less well studied. However, 
the literature does exhibit several quests for a semicoherent 
metric ifTTl [TSl l24l |251. About a decade ago, the first gen- 
eral discussion of the semicoherent metric for CW searches 
was given in lfT6l . along with numerical investigations. But 
only recently, the first fully analytical semicoherent metric has 
been found, leading to a significantly improved CW search 
method llT9l . This recent progress is based on a better un- 
derstanding of the global parameter-space correlations |26|, 
which were first examined in |27|. In turn, this insight pro- 
vides new coordinates on parameter space, enabling the ana- 
lytical calculation and study of the semicoherent metric. 

The present paper extends the recent work of |19| to 
greater generality and provides the essential technical basis 
for a parameter-space metric formalism for semicoherent CW 
searches. Additionally, complete and fully analytic semico- 
herent metric results are presented, which are ready to use for 
practical implementations of semicoherent searches, serving 
the earlier-mentioned quests of previous literature. The re- 
sults apply to broadband all-sky surveys (embedding directed 
searches with fixed sky position) for isolated CW sources. 

Section Inlbriefly reviews matched filtering for CW signals 
and the parameter-space metric formalism in general. Sec- 
tion |lll] elucidates how the metric is obtained for semicoher- 
ent searches. To evaluate and study the semicoherent metric, 
new coordinates on parameter space are defined in Sec. [TV 



In Sec.lv] the semicoherent metric is derived and investigated 
for the case of CW signals whose intrinsic frequency changes 
linearly with time (parametrized by one spindown parameter). 
In particular, the refinement factor is introduced as the ratio 
of the number fine-grid and coarse-grid templates, quantify- 
ing the additional parameter-space metric resolution due to 
combination of many segments. Section VI extends these re- 



sults to the case of CW signals whose intrinsic frequency can 
change quadratically with time (considering two spindown pa- 
rameters). Finally, a short conclusion follows in Sec.|Vn[ 



II. 



MATCHED FILTERING FOR CONTINUOUS 
GRAVITATIONAL-WAVE SIGNALS 



The detector output data time series is denoted by x{t) at 
detector time t. In the absence of any signal, the data con- 
tain only noise n{t), which is assumed to be a zero-mean, sta- 
tionary, and Gaussian random process. In case a signal h{t) 
is present, the noise is assumed to be additive, such that 
x{t) = n{t) + h{t). 

The dimensionless signal response function h{t) of an in- 
terferometric detector to a weak plane gravitational wave in 
the long-wavelength approximation is a linear combination of 
the form 



h{t) ^ F+{t)h+{t) + F^{t)h^{t), 



(1) 



where F+_x are called the beam-pattern functions, resulting 
in the amplitude modulations from Earth's spinning motion. 
They lie in the range — 1 < i^+,x < 1^ and depend on the 
orientation of the detector and source, and on the polarization 



angle ip of the waves. For explicit expressions the reader is 
referred to |13j. In the case of an isolated rapidly rotating 
neutron star with nonaxisymmetric deformations and negligi- 
ble proper motion (cf. f28','29 1), the waveforms corresponding 
to the plus (+) and cross (x) polarizations are 

h+{i)^ A+s\n^!(t), /ix(i) = v4x cos*(t), (2) 

where A+ and Ay, are the plus and cross polarization ampli- 
tude parameters, respectively, and '^{t) is given by 
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where $o is the initial phase, f^^^ = f denotes the frequency, 
and /('=>o) is the fcth frequency time derivative (also called 
"spindown"), evaluated at the SSB at time tg- The integer s > 
denotes the number of frequency time derivatives to be taken 
into account, therefore it holds f('^>^) — 0. The vector r{t) 
connects from the SSB to the detector, c is the speed of light, 
and ft is a constant unit vector pointing from the SSB to the 
location of the CW source. Thus, a point in phase parameter 
space p e T' is denoted by p = {f'^''\fij in respect of the 
reference time Iq. 

As first shown in |T3l, the phase ^(t) in Eq. ^ can be 
approximated without significant loss in signal-to-noise ratio 
(SNR) to good accuracy by 
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Consider a data segment spanning the interval 
[-T/2,T/2]. The J'-statistic is obtained ||T3] [SO) from 
the likelihood ratio A, which takes the form 



lnA={x\h)-^{h\h), 



where the scalar product has been defined as 

Hf)r{f) 
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(5) 



(6) 



with the Fourier transform indicated by"; the complex conju- 
gation denoted by *, and Sn defined as the one-sided noise 
spectral density. One may assume 5„ to be constant over the 
narrow bandwidth of the signal considered in this work. Then 
the scalar product of Eq. (J6]l is approximately given by 
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(7) 



Thus, the time average is introduced by the following nota- 
tion: 
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Using this notation, Eq. (J5]) is rewritten as 



lnA = 



2T 



{xh)--{e) 



(9) 



Replacing the amplitude parameters {yl-)_. ^x, ^, ^0} by 
their values which maximize In A, the so-called maximum 
likelihood (ML) estimators, defines the detection statistic F 
as 



J" = lnA 



ML- 



(10) 



Because the J^-statistic is maximized over the amplitude 
parameters, the remaining search space is just the phase- 
parameter space V. 



^sig(0 be described by known phase parameters denoted by 
the vector p^j^ and defined at t^, 

h,ig{t) = yli,sig sin $sig(<) + ^2,sig cos $sig(0 • (14) 

Since the signal parameters are known, one can construct a 
template which perfectly matches the signal. Assuming the 
noise n{t) to be stationary, Gaussian, zero-mean, and additive, 
one can show ll32l l33l that for a known signal (zero parame- 
ter offsets) the expectation value and variance of J-'*, respec- 
tively, are 



1 



E[T*] = l + -p'{0) 



'F" 



l + p'(0), 



(15) 



A. Coherent detection statistic for a simplified signal model 

Since the primary goal of this work is in relation to 
template-grid construction, a very useful approximation for 
this purpose is to replace the beam-pattern functions F+.x (i) 
by constant values IJ^Hll^S * 31 1. The phase of the CW signal is 
expected to change very rapidly at the terrestrial detector site 
over a time scale of typically less than ten seconds, whereas 
the amplitude of the signal varies with a period of one sidereal 
day. Coherent observation times of practical interest are typ- 
ically longer than one day, so that replacing the beam-pattern 
functions F+ x (i) with effective constant values is a good ap- 
proximation. In this case the signal model in Eq. (fT]| takes the 
form 

/i(i) = Ai sin$(<) + ^2 cos$(t), (11) 

where Ai 2 are defined to be the constant amplitude param- 
eters. In fSSl, the validity of this approximation is also in- 
vestigated using Monte Carlo simulations. It should be noted 
that the actual computation of the J^-statistic in a CW search 
will, of course, include the effects of amplitude modulation 
and involves precise calculation of the detector position with 
respect to the SSB using an accurate ephemeris model. This 
simplified signal is used here to facilitate the template-grid 
construction. 

The log likelihood of Eq. (J9|l for the simplified signal model 



( 11 1 is denoted by In A* and takes the form 



lnA* = 
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(12) 
By substituting the ML estimators for A12 in In A* of 



Eq. (12 1, it is straightforward to show |26, 27 1 that the sim- 



plified signal model (111 leads to the following coherent de- 
tection statistic J-* approximating F as 
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B. Perfect match of signal and template phase parameters 



Consider a signal h^igit) following the model of Eq. (Ill 
present in the data x{t). Let the signal's phase evolution 



where p(0) defines the optimal SNR, obtained as 
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The expression of Eq. ( [T6| may be further simplified using 
Eq. ^^ to yield 
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C. Mismatch of the signal and template phase parameters 

If the signal parameters are unknown in advance, one has to 
evaluate the detection statistic for a bank of templates. Let the 
template phase-parameter vector be p and the corresponding 
phase be $(t). In general, the template phase parameters will 
not exactly match the signal parameters. Thus, the parameter 
offsets are labeled by 



^P = Psig P 



(18) 



The resulting difference in phase A$(i) between the phase 
^sio{t) of the signal and the phase $(i) of a template is defined 

as 



A$(t) ^ ci>,,g(i) _ ci>(i) 



(19) 



In this case the expectation value and variance of J^*, respec- 
tively, are given by 

E[J^*] = l+^-p^Ap), a^.. =l + p2(Ap) (20) 

where the SNR p(Ap) here depends on the parameter off- 
sets Ap, such that 



AT 9 

p2(Ap) = — |(/.,„e-*)f 



(21) 



Further simphfication of Eq. (21 1 leads to 
T 
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The above relation shows that given an offset Ap between the 
signal and template parameters the squared SNR is reduced 
by I (e'^*) I . This gives rise to define a dimensionless "mis- 



match" M as 



M 
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The mismatch represents the fractional loss in the expected 
detection statistic due to the parameter Ap and thus provides 
a distance measure in the template parameter space. 



D. Metric on parameter space 

Taylor-expanding the mismatch A4 up to quadratic order in 
terms of the template parameter-space location offsets p at the 
signal location p,; yields 



x«^5,fcAp"Ap^ 

a.b 

defining a positive definite metric tensor g as 

gab^{da<^db<i>)~{da<^){db<i>) 



(24) 



(25) 



where a and b label the tensor indices, and the following no- 
tation has been employed: 



da^ 



(9p" 



(26) 



The expression of Eq. (25 i is often called the "phase metric" 
lfT5ll23ll29l . since it describes a distance measure on the phase 
parameter space V. 



ni. PARAMETER-SPACE METRIC FOR 
SEMICOHERENT SEARCHES 

In semicoherent CW searches the data are divided into N 
segments of duration T, where T is much smaller than one 
year. This allows to analyze each segment coherently, using a 
coarse grid of templates. Then the coherent detection statistic 
results from all segments are incoherently combined using a 
common ^ne grid of templates. This scheme is often called 
a semicoherent search strategy, offering the best overall sen- 
sitivity at limited computational resources fTS", Tl] when the 
fully coherent approach is infeasible. In preparation of calcu- 
lating the semicoherent metric, this section introduces some 
general notation. 



A. The coherent metric for a given segment 

Let the integer j = 1, ..., A^ label the jth segment, and let tj 
denote the time midpoint of segment j. The time average over 



the jth segment is defined by 
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The mismatch A4j in the jth segment is given by 
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where 
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Thus A4j can be approximated by 



■^.^E^afcAp'^Ap 



(27) 



(28) 
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(30) 



a,b 



where the components of coherent metric tensor for the jth 



segment g"' are obtained in analogy to Eq. (25 i as 



M 



gl^^ = {da^db^)y^^{da<^)y]{db^)y^ 



(31) 



When searching a subspace S-p of the phase parameter 
space V, the corresponding proper volume V is given by 



V = I dV = VdetgWdp. 



(32) 



The placement of signal templates to cover the search param- 
eter space is an instance of the sphere covering problem ll34l . 
Using a lattice of templates, the number of coarse-grid tem- 



plates Nt 
as 



is obtained from the coherent metric tensor g^^ 



PoV 



po I Vdct^bldp, 



(33) 



where the constant po describes the density of templates. The 
specific value of po depends on the desired maximum mis- 
match and the chosen type of lattice ifTSl [32l [34l . When using 
a random template bank instead of a lattice, then po can also 
depend on the desired coverage fraction 1.35,, 36J. 



B. The semicoherent metric for combining segments 

In the semicoherent search approach, coherent detection 
statistic results from the different segments are incoherently 
combined. To evaluate the metric for this case, the simplified 



coherent detection statistic T* (approximating T) of Eq. ( 13 1 
is again used here. Thus, J^* means the J^*-statistic value ob- 
tained in the jth segment. Recall that T* is the log likelihood 
function (maximized over the amplitude parameters). As the 
joint likelihood is the product, the joint log likelihood of all 
segments is the sum over j. Therefore, we define the semico- 
herent detection statistic J^ by 
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For the case of a known signal (zero parameter offsets) the 
expectation value and variance of T , respectively, are 
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(35) 



assuming identical noise spectral densities 5„ in every seg- 
ment. Hence, combining detection statistics from the N seg- 
ments reduces the variance by N. 

For nonzero offsets Ap between the template and signal 
parameters the resulting expectation value and variance of T , 
respectively, are obtained as 
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Thus, the mismatch A4, which measures the fractional loss 
in the expected semicoherent detection statistic T due to 
phase-parameter offsets Ap is obtained as 
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Since A4j is the mismatch in segment j according to Eq. ( (30| , 
Ai represents the average mismatch across the segments (cf. 
|fT6l ). Consequently, one may write 



X«^5,fcAp"Ap^ 



where the components of the semicoherent metric tensor g 
are obtained as the average of the individual-segment coherent 
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metric components g^^ from Eq. (Bib, 

1 ^ 
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(40) 



plates Ni is given by 



Thus, in analogy to Eq. (33 1 the number of fine-grid tem- 



AAfi- = Po / v/det|dp. (41) 



IV. NEW COORDINATES ON PARAMETER SPACE 



The analytic evaluation of the semicoherent metric compo- 
nents cjab from Eq. (|40]l is one of the central aspects of this 
work. This problem is approached by introducing new coor- 
dinates on the phase parameter space V, leading to a phase 
model which depends linearly on the coordinates. 



A. Linear phase model 

For coherent segment lengths T much smaller than one 
year, the orbital component rorb of the Earth's motion, 
r — Toib + ?^ipin, varies slowly during T and thus can be Taylor 
expanded around the segment's midpoint. Hence, by separat- 
ing the orbital and spinning components of the Earth's motion 
in the phase model $(t) a convenient reparameterization is 
obtained in which $(t) depends linearly on the new coordi- 
nates. For further details the reader is referred to Ref [26J . 
Thus the resulting phase model $(i) is obtained as 



$(t)=0o(to) + E^^'^(^o) 
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+ nx(io) cosil t + ny(fo) sin fit, 



(42) 



where 0o(io) is a constant independent of t, and the new fre- 
quency and frequency-derivative coordinates u'^^^to) as first 
derived in 1261 are 



v^'^\t^) = 2^ 
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(43) 



(39) and the new sky coordinates (as in lfT9l[3Tl ) are given by 



nx(to) = 27r/(to) te cosSd cosS cos [a — ao{to)] , (44a) 
ny{to) = 2nf{to) te cos Sd cos 6 sin[a — a_D(to)] • (44b) 

Thereby, £,{t) = foih{t)/c, with rorhit) denoting the 
vector from the Earth's barycenter to the SSB, and, 
Te = Re/c ~ 21 ms is the light travel time from the Earth's 
center to the detector, an (to), Sd are the detector position 
at time tg, and fi = 27r/(l sd) is the angular velocity of the 
Earth's spinning motion, which has a period of one sidereal 
day. 

Apart from an overall factor, the quantities i/^'"' have been 
referred to as the "global-correlation parameters" |fT9ll26l . As 
similarly done in earlier work 1231 [3 lil . note that the parame- 
ters i/'*''^ include a rescaling factor of (T/2)'^+^, such that they 
become dimensionless for convenience. 



The standard "physical" coordinates on V are the fre- 
quency and frequency derivatives /'■'^■'(to) at reference time 
to, and the unit vector n = (cos i5 cos a, cos S sin a, sin 6) on 
the two-sphere S'^, pointing from the SSB to the source. The 
angles a and 6 are right ascension and declination. 



B. Validity estimation 



The phase model of Eq. ( 42 1 is an approximation to the 
phase evolution described by Eq. (J4|i. The validity of this 
approximation depends on the coherent integration time T 



(duration of a given segment) and on the source frequencies 
searched, as was previously investigated in 1261 l?I |. The max- 
imum value of T as a function of the highest search frequency 
/ may be estimated by considering the first neglected term in 
the phase model becomes large enough to eventually lead to 
a significant mismatch. Given the phase model of Eq. (|42]i 
and searching s spindown parameters, the first neglected term 
A(j)s is given by 

A0,(t) ^ 2^ ^V^°rr fito)&+'Hto) ■ n (45) 
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The mismatch produced by the phase offset A(f>s follows from 
Eq. (28 I as 1 — [(c"^"^-")! . Figure 1 shows the values of T as 
a function of frequency / for which Acj)^ yields a mismatch 
of 30%. For instance, with s = 1, one should be able to use 
coherent segment durations T up to about 2 days for search 
frequencies / up to about IkHz, such that the mismatch due 
to the approximate phase model of Eq. (42 1 is still less than 
about 30%. These results also qualitatively agree with the 
earlier investigations reported in 1311 . 

However, most importantly, the new coordinates 
{t/('^\nx, riy} have significant advantages over the orig- 
inal coordinates {f^''\fi}. These new coordinates permit 
the first analytical solution for the metric of the incoherent 
combination step in hierarchical searches, yielding the "semi- 
coherent metric" to be presented in the following. In addition, 
in these new coordinates, the obtained metric is explicitly 
coordinate independent, a very convenient feature when it 
comes to practical aspects of conducting CW searches. 
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FIG. 1: Validity estimation of the approximate phase model Eq. l |42| l 
in terms of the number of spindown parameters s considered, for 
given values of coherent integration time T and the source frequen- 
cies searched. The dashed lines correspond to 30% mismatch for 
the one and two spindown case, respectively. For example, when in- 
cluding one spindown parameter (s = 1), the mismatch due to the 
approximate phase model should be less than 30% for coherent in- 
tegrations T up to about 2 days while searching frequencies up to 
about IkHz. 



V. METRIC EVALUATION FOR ONE SPINDOWN 
PARAMETER 

For all-sky surveys of prior unknown CW sources the 
search parameter space S-p is typically a four-dimensional 
subspace of V CHSl [IS, restricting to linear changes in 
frequency (one spindown parameter). In this case, using 
the new coordinates a point in V is labeled by the vector 
p = (i^, j>,T7,x, %) at a given reference time to. The phase 
model of Eq. (|42]) with s = 1 takes the form 
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where the coordinates f and v from Eq. ( 43 i are explicitly 
written as 
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y{h) - 2ti 



T 



f{to) + f{to)ato)-n + f^{to)-n 



T 



f , f{to) 



^ito)-n+fato)-n 



(48a) 



(48b) 



The coordinates n^ and Uy are as given by Eqs. (44i 



A. Coherent metric 



Using the coordinates {v,i>,n^,ny}, the components g]^l 
of the symmetric coherent metric tensor can be computed an- 



dimensionless quantity (p as 



alytically from Eq. pl\ . Thereby it is useful to define the 

(49) 



ip = n ■ 



T 



[j] 



Thus from Eq. (^3 111 the components g^^ are obtained as 
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"3' 




4 (t,-to\ 
3 [ T J ' 




4 16 ftj-to\ 
45 3 \ T J 




= - JiCv') sin (ntj) , 




= ji(v3) cos{ntj) , 



-4jiM ("^ 
9^-1^ -I hiv) siniQt,) 



sin (il t 



J/ 1 



+ 4jiM 



tj — to 

T 



cos i^tj) , 



(50a) 
(50b) 

(50c) 

(50d) 
(50e) 

(50f) 
(50g) 



1 1 



Jj] - i _ 



Ji] - _ 



jo{ip) cos{ip)~ ji{ip) sin (if) cos^ (Qtj) 



(50h) 



9nin, = - JiM sin(^) sm{ntj) cos {Qtj) , 



(50i) 



1 1 



fll^'l = - - 



jo(v5) cos(¥3)- ji(93) sin((^) sin^(fitj) 



(50j) 



where the spherical Bessel functions jn{x) [37] are defined 
by 



Ux) ^ i-xT 



1 d 



sin(a;) 



X dx / X 

The first few spherical Bessel functions are given by 

sin(x) 



Jo(a;) 

J2(a;) 
J3(a;) 



sin(x) cos(a;) 
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x^ 
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sin(a::) 3 cos(x) 



X 

sin(a::) 

X 



15 



cos(x) 



(51) 

(52a) 
(52b) 
(52c) 
(52d) 



Note that the components g^jj of the coherent metric tensor are 
explicitly independent of the coordinates. Therefore, the num- 
ber of coarse-grid templates J\f^°^'''^'^ as described by Eq. ( p3) l 
can be rewritten as 



Mr'"" = Po VdetgU] / dp , 



(53) 



where ^/detg^ has been taken outside the integration over 
the searched region of parameter space, since it is independent 
of the coordinates. Thus, ^/detg^ directly scales the num- 
ber of templates 7VJ'°'"'"'. The actual value of TVJ'"'"'"' depends 
on the parameter-space region S-p searched over To analyt- 
ically obtain realistic estimates for A/"™"^"^, one may assume 
the ranges 



n2 



<u< 



2 



~^~A ^ ^ ^ '^'a 



(54) 
(55) 



where Tmin = /// represents the "minimum spindown 
age" IIT3I to search for At fixed frequency / the parameter 
ranges of n^ and riy determine a two-dimensional disk Vf 
with radius of about 211 Jte- Thus, Eq. (53 1 yields 



Mr 
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PoVdetTgM 



drix driv 
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E rpS 



4 r,^ 
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(56) 
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FIG. 2: Dependency of the number of coarse-grid templates _/V™"™ 
on the coherent integration time T (segment length) for the one- 
spindown case. The solid curve shows T^ ydetgLil, since 
y^coarse ^ rp3 ^detgbl. The dashed curve represents rV\/l35. 



The determinant of the coherent metric tensor det 17 [^1 is 
obtained as 



det g[^l 



1 
135 



[l-6ji2(^)- jo((/7) cos((p)] 



X [l-10J2^((/?)- jiM sin((^) 

- jo(¥?) cos((y9)] . (57) 

Note that det gl-'l is also explicitly independent of the refer- 
ence time to (as well as of tj), and solely depends upon ip. 
This can be understood from the following reasoning. For the 
parameters /^^^^ a change of reference time corresponds to a 
linear transformation |28, 29|, whose determinant is 1. Simi- 
larly, for the parameters n^ and riy a change of reference time 
can be described by a rotation, whose determinant is also 1. 
Therefore this explains why det gl^l is independent of the ref- 
erence time. 

To examine the scaling of ^'=°'»rse ^j^j^ j-j^g coherent inte- 
gration time T, Fig. pi shows T^ -\/det gbl, which is a mea- 
sure for Af™''"^^'^, as a function of T. For increasing values 
of T, detgl^l converges to 1/135, since the metric tensor 
components related to the Earth's spinning motion, n^ and 
riy, become approximately constant. This behavior can be 
understood from the Rayleigh criterion. An estimate of the 
diffraction-limited resolution is described by the ratio of the 
wavelength and the effective "baseline" ["^Jl. The maximum 
baseline in terms of the Earth's spinning motion is the Earth's 
diameter, which is first reached already after half a day of inte- 
gration. Therefore, very little metric resolution is gained after 
integration times T of about a day, as long as the Earth's or- 
bital motion can still be well modeled by a Taylor expansion. 

Finally, it should be mentioned that a convenient choice 
of to in favor of a compact notation is to = tj = 0. To indicate 
this choice has been made, the resulting coherent metric ten- 
sor is denoted by gf-'^"'. The components g^^^ are obtained 



from Eqs. (50 1 as follows: 



using the following notations to simplify the expressions: 





1 




= 0, 




4 
~45' 




= 0, 




= JiM, 




= -3 J2('^) 


UVUy 


= 0, 




= 2^2^"^' 


i/rixiiv 


= 0, 



(58a) 

(58b) 

(58c) 

(58d) 
(58e) 

(58f) 

(58g) 

jo(¥') cos((/?) - ji((^) sin(v?) , (58h) 



gliJl = 2 " 2^°''^'* cos((/?) 



(58i) 
(58j) 



B. Semlcoherent metric 



Given the coherent metric tensor components g^l in 
Eqs. dSOll, explicit expressions for the components cjab of the 



semicoherent metric tensor g are obtained via Eq. (40 1 as 
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_ _ 4 16 
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9>^ny =Ji(v)Mo°% 



?.«, = -3 j2MAir+4jiMMr 
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(59f) 
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ffnx«x = 2 " 2-'°'''^'' *^°'^ *'^'* ~ ^^^"^^ ™ '^'^'' ''™' ' (59h) 
5nxn. = -JiMsin(^)Cr™% (59i) 

5"y«y = 2 " 2-'°'-'^'* '^°'^ *'^'' " Jl'^'^) '^^^ ('^^ ''2™ ' "^^^J^ 
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tj — to 
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N 
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1 ■^^ f tj — to 

N ^ 

1 ^ 

-T 



N 



tj — to 



sin {ntj) , 



cos(r2ij) , 



(60) 
(61) 
(62) 



w JV 



J=l 



AT 



J = l 



Cr°^^^EBin(17i,)cos(J7i,), 



(63) 
(64) 



j=i 



where m can be zero or a positive integer number. 

The components (jab of the semicoherent metric tensor are 
also explicitly independent of the coordinates. Therefore, the 
number of fine-grid templates A/^*^™ as described by Eq. (41 1 
is rewritten as 



AA«"'= = po Td^ / dp 



-K" 



I detg 
det gbl 



(65) 



Considering the distribution of the segments' mid- 
points {tj}, the quantities /i.,„ may be interpreted as the mth 
moment of this distribution. Thus a very convenient choice of 
reference time to for the semicoherent metric is given by the 
time average of all segment's midpoints {tj}. 



1 ^ 






(66) 



With this choice of to the quantities /!.,„ become the mth cen- 
tral moments (denoted by /!„) of the distribution of the seg- 
ment midpoints {tj}. 

In addition, for a large number of segments N, where the 
data set spans over many cycles of the Earth's spinning motion 
with the one-day period 27r/rj, the following approximations 
may also be used: 



TV 



^^sin(nt,)«0, 
1 ^ 1 

Y^sin'int,)^- 



1 ^ 
-Y,cos{ntj)^0, (67a) 



N 



3 = 1 



N 



l^cos2(f^i,)«i, 



3 = i 



— Y^ sin {n tj) cos (rj tj) « , 



(67b) 
(67c) 



j=i 



as well as 

N 



N 



t. 



(67d) 
Hence, given ^™ « 0, ^^°' w 0, Atf' ~ 0, /if' w 0, 
CT w 1/2, C™' « 1/2, and Cf™™' « 0, the semicoherent met- 



ric components of Eqs. (59i along with the tg choice of 



Eq. ( 66 1 take the following diagonal form: 
/I 
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^ 45 






fh 






Rjv) 
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2 



(68) 



where R{tf) is defined as 

^('^) ^ 1 - J ('/') cos ((p) - j 1 ((/?) sin {if) . (69) 

The determinant of the above semicoherent metric tensor gab 



from Eq. (68 1 is obtained as 
det^ 



1 
135 



9 



M2 



R\ip). 



(70) 



Finally, if T is an integer multiple q of one sidereal day, 
T = ^g, such that ip — irq and R{Tiq) = 1, the metric ten- 
sor g from Eqs. (68 1 simplifies to 



(\ 







o\ 



" 45 ' 3 





i^A2 















(71) 



M 



and the corresponding determinant is simply given by 

1 4 . 



detg 



135 9 



M2. 



(72) 



C. Parameter-space resolution refinement 



The results for the semicoherent metric tensor shown in 



Eqs. (68 I feature an important property: gij^ represents the 
only component of the semicoherent metric tensor g which 
significantly changes with an increased number of data seg- 
ments N , not converging to some constant value. In this re- 
spect, increasing N means that the number of fine-grid tem- 
plates needs to be increased in only one dimension compared 
to a given coarse grid. 

To describe the refinement quantitatively, we use Eq. (65 i 
to introduce the refinement factor, denoted by 7, defining the 
ratio of the fine and coarse template-grid points: 



_yy-coai 



/ detg 
detgbl ' 



(73) 



where in the last step Eqs. (53 1 and (65 1 have been used. 
In what follows, refinement factor 71 (72) with the sub- 
script 1 (2) is used to indicate the one-spindown (two- 
spindown) case, respectively. 



Thus, by means of Eqs. (57i and (70i the refinement fac 



tor 71 for the one-spindown case is explicitly given by 



71 = V1 + 6OA2 Qiv) , 

where the function Q{<f) is defined as 

Q{^) = Rip) [I - 6 ji2(^) - joM cos M]"'/' 
X [1-10J2^(V')- JiM sin(^) 
- jo((^) cos((^)]^^ ^ 



(74) 



(75) 



Note that from Eq. (74i it is obvious that 71 scales linearly 



with the number of data segments N, since only JI2 de- 
pends (quadratically) on TV. Hence, the enhanced parameter- 
space resolution resulting from the incoherent combination in- 
creases approximately as ex A^ solely due to the first spindown 
parameter This is related to the fact that the number of possi- 
ble (linear) spindown tracks in frequency across all segments, 
of course, grows linearly with N, too. 



D. Illustrative example 

As for a simple example, one may consider the case 
where N segments are uniformly distributed in such a way 
that there are no gaps between neighboring segments, so 
that tj = [j - (N + l)/2] T. It should be pointed out that 
this special case had been assumed a priori in the previous 
work of lfT6l . Thus, for this particular instance the numeri- 
cal findings of fTSl can be compared to the analytic results 
found here. Denote the time span of the entire data set as 
2data = NT. Thus, in the present example the choice of 



Eq. ( 66 1 yields to = and 



N^ 



/"2 



12 



(76) 



For this case, 71 is shown in Fig. l3] for different values of 
T and A^. For increasing values of T, 71 rapidly converges 
to some upper-limit value for fixed N. This maximum con- 
stant value can be arrived at analytically in the following way. 
In this regime of T, det gl^' from Eq. (57 1 can be approxi- 
mated by 1/135 and det g by Eq. (72 1, leading to Q{(p) ~ 1 
in Eq. ( [74| and thus 



71 « Vl + 60/i2 = \/5N^ ~ 4 . 



(77) 



where Eq. (76 1 has been used in the latter step. This result 
agrees well with the corresponding refinement factor that can 
be computed from the numerical template counting formulas 
given in |,16l. 
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FIG. 3: Refinement factor 71 for the one-spindown case. In this plot 
the color-coded contours show logiQ(7i) as a function of the coher- 
ent segment length T and the number of segments A'^. It has been 
assumed that there are no gaps between neighboring segments, and 
so Eq. J76| has been used. The dashed lines mark locations of data 
sets with total time span (Tdata = NT) of one, two, and three years. 



VI. METRIC EVALUATION FOR TWO SPINDOWN 
PARAMETERS 



Searching for CW signals from potentially younger objects, 
spinning down faster, might eventually require one to also in- 
clude a second spindown parameter (see 1*381 for an exam- 
ple), taking also into account quadratic changes in frequency 
with time. In this section, the additional semicoherent metric 
components are computed, which arise when a second spin- 
down parameter is included in the search. Thus, the coordi- 
nates {h', j>, t>, fix, Uy} at time Iq are used to label a point in 
this five-dimensional parameter space. The phase model of 



Eq. ( 42 1 with s = 2 is given by 






T 



212 



8 (t — t )^ 
i^ih) 7p^ ^' n-x(io) cosili + ny(to) sinfJi. 



T3 



(78) 



Based on Eq. (43 1 the coordinates v, v, and v are explicitly 
written as 



i/(io) - 27r 
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/(io) + /(io)C>o)-^ 



Ai^)iih)-n 



(79a) 
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f{to) , f{to) 



e(to) • n 



+ f{to)^ito)-n+^ato)-n 



(79b) 



i>(io) = 27r 



T 



f , /fa) 
6 6 



e(io)-n 



+ ^'iito)-n+^i{to)-n 



(79c) 



The coordinates n^ and Uy are as introduced in Eqs. (44i 



A. Coherent metric 



Including a second spindown parameter i> in the phase 



model of Eq. ( 42 1 yields the following additional components 
for the coherent metric tensor gl^': 

/ \ 2 

„bl _ 1 I /I / *J ^ ^0 
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tj — to 
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2 

sin {Qtj) 
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(80a) 
(80b) 
(80c) 



(80d) 



^) cos(f]t,) 



(80e) 



The components gj^jj of the coherent metric tensor are ex- 
plicitly independent of the coordinates as mentioned earlier 
Therefore, the number of coarse-grid templates _/VJ=°^"e j^ ^j^^ 
computed as presented by Eq. ( [53] l. To analytically estimate 
the actual value of _/VJ^°''rse Jqj. jjjg two-spindown case, in ad- 



dition to Eqs. ( |54| ) and ( |55] l the following ranges of i) are as- 

TV 



sumed to be searched 



12r 2 

J^ -^ ' min 



< D < 



-•-^ 'mill 



(81) 



10 



Thus, evaluation of Eq. (53 1 yields in this case 



povdetgW 



dv 



^T/„ 



driv driv 



tT^J 



dv 



Av 



poVdetgW 



6^2 



TT r 



T' (/ix - /^.„) ■ (82) 



The determinant of the coherent metric tensor gl-'l for the 
two-spindown case is obtained accordingly as 

4 



det qf-'l = 



[l + 375ji^((/.)J2'M + 189ji"M 



23625 

-252ji3(^)J3(^) + 84ji2(^)J32(^) 

+ 42ji(^) J3((/7) -42jo((^) jiM J3(</?) cos((/7) 

+ 75ji3(<^)sin(^)/2-14J32(^)-10J22M 

+ 14joMJ3^(<y5) cos{Lp)-m2i^{if) 

+ 10 jo('^) J2^(V5) cos((^) - 2 jo((p) cos((p) 

+ 69joMji^(<p) cos((/3) + jo2((/3) cos2(^) 

- ji((^) sin(^) + jo((y3) ji((^) sin((y3) 003(93)] . (83) 

As mentioned earlier, detg'-'l only depends upon cp. To in- 
vestigate the scaling of _/V]^°a"<= with the coherent integra- 
tion time r. Fig. MJ shows T^ydctgl^l, being a measure 
for A^"""'', versus T. For increasing values of T, detgl-'l 
converges to the constant value of 4/23625, since the met- 
ric tensor components related to the Earth's spinning motion, 
n^ and riy, become approximately constant in this regime as 



explained earlier in Sec. V A 



Finally, it should be pointed out that choosing Iq = tj = 
simplifies the expressions of Eqs. ([80|, yielding 






■ jilv?)- E J3('y3) 



(84a) 
(84b) 
(84c) 
(84d) 
(84e) 



B. Semlcoherent metric 

The extra components of the semicoherent metric tensor g 
I usii 

9vv = ^ +'^^J'2 



are then obtained via Eq. (40 1 using Eqs. (80 1 as 
1 



duu 



Ml + 16Ai3, 



9vv=Tj+^ fJ-2 + 48 /X4 , 
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(85b) 
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(85d) 
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With the approximations of Eqs. (67 1 used earlier, /ig 

0, /if 



f^r 



0, 

0, /i™' w 0, and the to choice of Eq. (|66]) the 



semicoherent metric tensor components in Eqs. (85 1 take the 
following form 



gi^u 


1 

^5 


goi) 


«0, 


goi} 


1 
"7 


gi>n^ 


« — 


QOUy 


«6j 



4^2, 



) /i2 + 48 /i4 , 



6jl(<^)A2°' ~ -6ji((/?)A2, 

6jiMAr«6jiMA2. 



(86a) 

(86b) 

(86c) 

(86d) 
(86e) 



Thus, the full five-dimensional semicoherent metric tensor g 
is obtained as 
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FIG. 4: Dependency of the number of coarse-grid templates _/VJ^°''"'= 
on the coherent integration time T (segment length) for the two- 
spindown case. The solid curve shows T^ \/detgbl, since 
y^coarse ^ y6 y'jet c/M . The dashed curve shows rV4/23625. 
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r 



where i?((p) is given by Eq. (69 1. The corresponding determi- 
nant is obtained as 



deig • 



l + 60/i2 
23625 



4i?2(^) 



X 140[l5(/i4-A2)-45 



.Ji%^A| 



M2 



(88) 



When T is an integer muhiple q of one sidereal day, 
T —'^q, such that ip — irq and R{nq) = 1, the determinant 
of the semicoherent metric tensor takes the form 



detg ; 



l + 60/i2 
675 



16[15(A4-A2) 



45 Ai 
■n'^q^ 



A*2 



(89) 



C. Parameter-space resolution refinement 



The refinement factor for the two-spindown parameter 



case 72 has been defined through Eq. ( 73 1 as 



72 



/ detg 
det gbl 



(90) 



Here g^^^ and g denote the coherent and semicoherent met- 
ric tensors, respectively, for the two-spindown parameter case. 
Substituting det g'-'l by Eq. ( 83 i and det g by Eq. ( 88 i yields 



72 = 2^35 yr 
15(A4 - 



- 60 A2 
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45 



Ji^(y)Ai 
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1/2 



u{^).. 



(91) 



where the explicit expression for C/((p) can be deduced from 
Eq. ( 83 1 (suppressed for brevity here). 



From Eq. ( 9 1 1 the scaling of 72 at leading order in iV is ob- 
tained as 72 oc N^, using A4 oc TV' and A2 oc N'^. This cubic 
scaling with N is solely due to the first and second spindown 
parameters, resulting from the product of possible linear and 
quadratic spindown tracks in frequency across the segments, 
which obviously grows as N^. 

Note that in general, for s spindown parameters, the ex- 
pected scaling of the refinement factor 7^ with the number of 
data segments N is given by 



7, oc 7V-(^+l)/2 . 



(92) 



This scaling with TV is robust and agrees with what can be 
deduced from the numerical findings of Ref . Ill 61 . 



D. Illustrative example 

To further examine 72, it is instructive to consider again 
the example data set presented in Sec. |VD| Thus, given 
tj = [i - (iV + l)/2] T and t^ = 0, the third central moment 
vanishes. As — 0> and the fourth central moment A4 is ob- 
tained as 
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+ 45 



(93) 



(94) 



Using Eq. ( 93 1 along with Eq. ( 76 1 the refinement factor 72 
can be computed as a function of T and A^, as illustrated 
in Fig. I5] For increasing values of T, 72 rapidly converges 
to some constant value for fixed N. In this case, detgl-'l 
of Eq. ( [83] l is well approximated by 4/23625 and det g of 
Eq. ^ by 
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FIG. 5: Refinement factor 72 for the two-spindown case. In this 
plot the color-coded contours show log 5^0(72) as a function of the 
coherent segment length T and the number of segments A'^. It has 
been assumed that there are no gaps between neighboring segments. 
The dashed lines mark locations of data sets with total time span 
(Tdata = NT) of one, two, and three years. 
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Hence, 72 in this case is described by 

72«y^2ioo(i + 60A2)(/i4-Ai) 

w 71 \/(35iV4_i75 7v2 + i40)/3 . (96) 



where Eqs. ( [76] l and ( [94] i have been used. Thus, the antici- 
pated scaUng at leading order in N is recovered: 72 oc N^. 



VII. CONCLUSION 

A formaHsm has been presented for the parameter-space 
metric of semicoherent CW searches, where the data are di- 
vided into segments that are coherently analyzed and subse- 
quently combined incoherently. By using new coordinates on 
parameter space, the first fully analytical semicoherent metric 
for broadband all-sky CW surveys has been derived. Addi- 
tionally, in the new coordinates the components of both the co- 
herent and the semicoherent metric tensor are constant, being 
explicitly independent of the coordinates. This entails great 
convenience regarding practical aspects of semicoherent CW 
searches. 

Explicit analytic expressions of the semicoherent metric 
tensor components have been obtained for two typical search 
parameter spaces of current practical interest. First, the one- 
spindown case has been considered, restricting to linear fre- 
quency drifts with time as done in many current all-sky 
searches. Second, the semicoherent metric also has been cal- 
culated and examined for the two-spindown case, where addi- 
tionally quadratic changes in frequency are taken into account. 

Analytic signal template counting formulas have been pro- 
vided for the coherent stage (coarse grid) as well as for the 
incoherent combination step (fine grid). In this respect, a 



useful quantity, called the refinement factor, has been intro- 
duced as the ratio of the number of fine-grid and coarse-grid 
templates. Thus, the refinement factor describes (coordinate 
independently) the additional parameter-space metric resolu- 
tion gained from the combination of segments. Moreover, the 
scaling of the refinement factor with the number of segments 
has been found to be predominantly determined solely by the 
spindown parameters. 

The present results also embed the case of directed semico- 
herent searches, where the sky position is known and hence is 
not a search parameter Thereby, the search parameter space 
consists only of frequency and spindown parameters. The cor- 
responding semicoherent metric tensor components are iden- 
tical to ones derived here. The resulting scalings of the re- 
finement factor with the number of segments also hold, since 
governed by the spindown parameters as described above. 

The formalism presented in this paper assumes segments 
that are very short compared to one year. The Earth's orbital 
motion then can be described by a low-order Taylor expan- 
sion, and therefore be modeled by changes in frequency and 
frequency derivatives. In the future, however, increased com- 
puting power might allow one to use segments substantially 
longer than a few days and hence issues remain to be explored 
in such circumstances. 



VIII. ACKNOWLEDGMENTS 

I am indebted to Bruce Allen, Reinhard Prix, and Chris 
Messenger for very useful discussions and helpful comments 
on this work. I also thank Miroslav Shaltev for corrections of 
the manuscript. I gratefully acknowledge the support of the 
Max-Planck-Society. This document has been assigned LIGO 
Document Number LIGO-P1000006-v4. 



[1] B. J. Owen, L. Lindblom, C. Cutler, B. F. Schutz, A. Vecchio, 

and N. Andersson, Phys. Rev. D 58, 084020 (1998). 
[2] G. Ushomirsky, C. Cufler, and L. Bildsten, Mon. Not. Roy. As- 

tron. Soc. 319, 902 (2000). 
[3] C. Cutler, Phys. Rev. D 66, 084025 (2002). 
[4] D. I. Jones and N. Andersson, Mon. Not. Roy. Astron. Soc. 331, 

203 (2002). 
[5] B. J. Owen, Phys. Rev. Lett. 95, 211101 (2005). 
[6] B. Abbott et al. (The LIGO Scientific Collaboration), Rep. 

Prog. Phys 72, 076901 (2009). 
[7] B. Abbott et al. (The LIGO Scientific Collaboration), Phys. Rev. 

D 77, 022001 (2008). 
[8] B. Abbott et al. (The LIGO Scientific Collaboration), Phys. Rev. 

D 79, 022001 (2009). 
[9] B. Abbott et al. (LIGO Scientific Collaboration), Phys. Rev. D 

80, 042003 (2009). 
[10] http://einstein.phys.uwm.edu/. 
[II] B. Abbott et al. (The LIGO Scientific Collaboration), The As- 

trophysical Journal Letters 683, L45 (2008). 
[12] B. Abbott et al. (The LIGO Scientific Collaboration), Phys. Rev. 

Lett. 102, 111102(2009). 
[13] P Jaranowski, A. Krolak, and B. F. Schutz, Phys. Rev. D 58, 



063001 (1998). 
[14] C. Cutler and B. F Schutz, Phys. Rev. D 72, 063006 (2005). 
[15] P R. Brady, T. Creighton, C. Cutler, and B. F. Schutz, Phys. 

Rev. 57,2101(1998). 
[16] P R. Brady and T. Creighton, Phys. Rev. D 61, 082001 (2000). 
[17] C. Cutler, I. Gholami, and B. Krishnan, Phys. Rev. D 72, 

042004 (2005). 
[18] B. Krishnan, A. M. Sintes, M. A. Papa, B. F Schutz, S. Frasca, 

and C. Palomba, Phys. Rev. D 70, 082001 (2004). 
[19] H. J. Pletsch and B. Allen, Phys. Rev. Lett. 103, 181102 (2009). 
[20] R. Balasubramanian, B. S. Sathyaprakash, and S. V. Dhurand- 

har, Phys. Rev. D 53, 3033 (1996). 
[21] B. J. Owen, Phys. Rev. D 53, 6749 (1996). 
[22] I. Jones, B. J. Owen and D. Whitbeck, LIGO Document 

T0900500-V1, 2005, available at https://dcc.ligo.org/. 
[23] R. Prix, Phys. Rev. D 75, 023004 (2007). 
[24] S. Dhurandhar, B. Krishnan, H. Mukhopadhyay, and J. T. Whe- 

lan, Phys. Rev. D 77, 082001 (2008). 
[25] A. L. Watts, B. Krishnan, L. Bildsten, and B. F Schutz, Mon. 

Not. Roy Astron. Soc. 389, 839 (2008). 
[26] H. J. Pletsch, Phys. Rev. D 78, 102005 (2008). 
[27] R. Prix and Y. Itoh, Class. Quant. Grav. 22, S1003 (2005). 



13 



[28] P. Jaranowski and A. Krolak, Phys. Rev. D 59, 063003 (1999). [34] R. Prix, Class. Quant. Grav. 24, S481 (2007). 

[29] P Jaranowski and A. Krolak, Phys. Rev. D 61, 062001 (2000). [35] C. Messenger, R. Prix, and M. A. Papa, Phys. Rev. D 79, 

[30] P Jaranowski and A. Krolak, Living Rev. Relativity 8, 3 (2005), 104017 (2009). 

http://www.livingreviews.org/lrr-2005-3/. [36] G. M. Manca and M. Vallisneri, Phys. Rev. D 81, 024004 
[31] P. Astone, K. M. Borkowski, P. Jaranowski, and A. Krolak, (2010). 

Phys. Rev. D 65, 042003 (2002). [37] M. Abramowitz and I. A. Stegun, Handbook of Mathematical 
[32] P. Astone, K. M. Borkowski, P. Jaranowski, M. Pietka, and Functions with Formulas, Graphs, and Mathematical Tables 

A. Krolak, Phys. Rev. D 82, 022005 (2010). (Dover, New York, 1964). 

[33] P. Jaranowski and A. Krolak, Analysis of gravitational-wave [38] K. Wette etal.. Class. Quant. Grav. 25, 235011 (2008). 

data (Cambridge University Press, Cambridge, 2009). 



14 



